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Abstract 

We review free fermion, melting crystal and matrix model representations of wall- 
crossing phenomena on local, toric Calabi-Yau manifolds. We consider both unrefined 
and refined BPS counting of closed BPS states involving D2 and DO-branes bound to a 
D6-brane, as well as open BPS states involving open D2-branes ending on an additional 
D4-brane. Appropriate limit of these constructions provides, among the others, matrix 
model representation of refined and unrefined topological string amplitudes. 
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1 Introduction 



This review is devoted to some aspects of counting of BPS states in a system of Dp- 
branes, with even p, in type IIA string compactifications. The problems of BPS counting 
span a vast area of research in supersymmetric gauge and string theories. Their important 
feature is a special, non-constant character of BPS multiplicities: their values depend 
on various moduli and jump discontinuously along some special loci in the corresponding 
moduli space, so called walls of marginal stability. The pattern of these jumps follows wall- 
crossing formulas, found from physical perspective by Denef and Moore [T] and, in more 
general context, formulated mathematically by Kontsevich and Soibelman P]. The regions 
of the moduli space in between walls of marginal stability, in which BPS multiplicities are 
(locally) constant, are called chambers. 

The BPS states we are interested in, and which we will refer to as closed BPS states, 
arise as bound states of a single D6-brane with arbitrary number of DO and D2-branes 
wrapping cycles of a toric Calabi-Yau space. More generally, we will also consider open BPS 
states, which arise when an additional D4-brane spans a lagrangian submanifold inside the 
Calabi-Yau space and supports open D2-branes attached to it. The closed and open BPS 
states give rise, respectively, to single-particle states in the effective four- dimensional and 
two-dimensional theory (in remaining, space-time filling directions of, respectively, D6 and 
D4-branes). In this context the character of BPS multiplicities can be understood in much 
detail, and it relates to other interesting exactly solvable models: free fermions, crystal, and 
matrix models. In brief, these connections arise as follows. Firstly, BPS states we consider 
turn out to be in one-to-one correspondence with configurations of certain statistical models 
of melting crystals. The structure of these crystals depends on geometry of the underlying 
Calabi-Yau space, as well as on the chamber one is considering. In consequence BPS 
counting functions, upon appropriate identification of parameters, coincide with generating 
functions of melting crystals. It turns out that the structure of these crystals can be given 
a free fermion representation. Furthermore, once such free fermion formulation is known, 
it can also be represented in terms of matrix models. Connection with vast theory of 
matrix models has many interesting mathematical and physical consequences and allows 
to shed new light on wall-crossing phenomena. The aim of this review is to explain these 
connections. 

The BPS generating functions which we consider are intimately related to topological 
string amplitudes on corresponding Calabi-Yau spaces. This relation is most transparent 
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in the physical derivation discussed in section [21 which rehes on hfting the D-brane system 
to M-theory. The M-theory viewpoint makes contact with original formulation of closed 
topological strings by Gopakumar and Vafa |3lll], and open topological strings by Ooguri 
and Vafa 0. In particular, in one specific, so called non- commutative chamber, the BPS 
generating function is given as the modulus square of the topological string partition 
function. In all other chamber BPS generating functions can be uniquely determined from 
that non-commutative result. There is also another special, so called commutative chamber, 
in which BPS generating function coincides (up to the factor of MacMahon function) 
with the topological string partition function. For toric manifolds which we consider, 
such topological string amplitudes can be constructed, among the other, by means of the 
powerful topological vertex formalism [6]. Relation to crystal models was in fact first 
understood in this topological string chamber jZl [H [^. One advantage of the formalism 
presented in this review is the fact that it allows to construct matrix model representation of 
all these generating functions (so, in particular, matrix model representation of topological 
string amplitudes). 

In more detail, we will consider generating functions of D2 and DO-branes bound to a 
single D6-brane of the following form form 

ZBPsiqs,Q) = J2^ia,/3)q':Qf', (1) 

where a G Z is DO-brane charge, and /3 G H2{X,'L) is D2-brane charge. Multiplicities 
VL{a, P) jump when central charges (which itself are functions of Kahler moduli) of building 
blocks of a bound state align, and therefore these generating functions are locally constant 
functions of Kahler moduli. Along the walls of marginal stability the degeneracies Q{a, (5) 
change and indeed obey wall-crossing formulas of lU |2] mentioned above. 

If there is an additional D4-brane which spans a lagrangian submanifold inside the 
Calabi-Yau space, in addition to the above closed BPS states, one can consider also open 
BPS states of D2-branes with boundaries ending on a one-cycle 7 on this D4-brane. In this 
case the BPS states arise on the remaining two-dimensional world-volume of the D4-brane. 
The holonomy of the gauge field along 7 provides another generating parameter z, so that 
open BPS generating functions take form 

ziTsiQs, Q) = Y1 ^("' (3)^:q^^'- (2) 

As we will show, generating functions of such open BPS states can be identified with 
integrands of matrix models mentioned above. 
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One more important aspect of BPS counting is referred to as refinement, and amounts 
to refining BPS counting by introducing one more parameter, customarily denoted /3. The 
refinement can be introduced from several perspectives which give rise to identical results, 
however their fundamental common origin is still not fully understood. We will introduce 
refinement by distinct counting of states with different SU{2) spins inside spacetime 5*0(4) 
rotation group in the generating function ([T]). In |1U] it was argued that this physical 
viewpoint should agree with the mathematical counterpart of motivic deformation [2], and 
also a refined version of a crystal model was constructed. Another notion of refinement 
arises in Nekrasov partition functions, which are defined in a non-trivial gravitational (so 
called Q-) background parametrized by two parameters ei and €2 [Uj. Nekrasov partition 
functions can also be defined for five-dimensional gauge theories and then they agree with 
topological string amplitudes. In particular, the formalism of the topological vertex [6j 
has also been extended to the refined context in [12j, and shown to reproduce relevant 
Nekrasov partition functions. Also BPS generating functions, in the limit of commutative 
chamber, are known to reproduce refined topological string amplitudes with /3 = — ei/e2 
|13] . However the worldsheet definition of refined topological string amplitudes is not fully 
understood. 

As an exemplary and, hopefully, inspiring application of the entire formalism presented 
in this review, in the final section [6] we derive matrix model representation of the refined 
topological string partition function for the conifold. The refined matrix model which we 
find has a standard measure, however its potential is deformed by /3-dependent terms. It is 
obtained by constructing appropriate refined crystal model and free fermion representation, 
and subsequently reformulating this representation in matrix model form. Finally, taking 
the limit of the commutative chamber, we obtain matrix model representation of the refined 
topological string amplitude. Even though we demonstrate this result in the conifold case, 
with some technical effort it can be generalized to other toric manifolds which we consider]^ 

Short literature guide 

The literature on the topics presented in this paper is extensive and still growing, and 
we unavoidably mention just a fraction of important developments. The relation between 

^As we recall in section [HI refined topological string amplitudes were also postulated to be reproduced 
by another type of matrix models, so-called /3-deformed ones (whose Vandermonde measure is deformed by 
raising it to power /3) ; however explicit computations showed that this cannot be the correct representation 
of refined amplitudes. 
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Donaldson-Thomas invariants for the non-commutative chamber of the conifold was first 
found by Szendroi [13]. It was generahzed to orbifolds of C^, and related to free fermion 
formalism, by Bryan and Young [15]. The relation to free fermions and crystals was 
extended to a large class of toric manifolds without compact four-cycles |16[ IT7] . These 
developments were accompanied by other mathematical works [TS| [T^ . 

In parallel to the above mentioned mathematical activity, wall-crossing phenomena for 
local Calabi-Yau manifolds were analyzed from physical viewpoint. The analysis of non- 
trivial BPS counting for the conifold was described by Jafferis and Moore in [20]. This 
and more general cases were related to quivers and crystal models in f2T\ 12^ . Derivation 
of BPS degeneracies from M-theory viewpoint and relation to closed topological strings 
were discussed in [23J, and generalized to open BPS counting in [211 |25l [261 EZ]. Relations 
to matrix models, discussed for plane partitions with some other motivation in [28|, were 
extended to other crystal models relevant for BPS counting in [29j, and also in [30j. Sub- 
sequently it was related to open BPS counting in [27j. Refined BPS counting was related 
to crystal models in [TOl [13], and corresponding matrix models were constructed in [31\ . 

Let us also mention some other, related works devoted to crystals and free fermions. 
The fermionic construction of MacMahon function for was originally presented in [7], 
and its relation to open topological strings and more complicated Calabi-Yau manifolds 
were discussed in [32 1 l33 l IM ] . Newer ideas, analyzing more complicated systems involving 
D4-branes, were presented in [35t l36] . More expository presentations of various aspects 
described here can be found in ^S7\ [55] . A general introduction to mathematical and 
physical aspects of mirror symmetry can be found in [39j. 

Plan 

The plan of this paper is as follows. In section [2] we introduce BPS generating functions 
and present one possible derivation of their form, which relies on the M-theory interpre- 
tation of a D-brane system, following [231 1211 [23 [2Z]. In section [3] we provide a little 
mathematical background and introduce notation pertaining to toric Calabi-Yau mani- 
folds, free fermion formalism, and matrix models. In section [H we introduce fermionic 
formalism for BPS generating functions and present corresponding crystal models, build- 
ing on earlier ideas of [3 [15] and following [^16j. In section [5] we reformulate the problem of 
closed BPS counting in terms of matrix models and relate it to open BPS counting [27 1 [2U ] . 
In section [6] we refine our analysis, present refined BPS generating functions and crystals 
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[10] , and construct corresponding refined matrix models [31] . 

2 BPS generating functions 

In this section we introduce generating functions of BPS states of D-branes in toric 
Calabi-Yau manifolds. Our task in the rest of this paper is to provide interpretation of 
these generating functions in terms of free fermions, melting crystals and matrix models. 
These generating functions can be derived using wall-crossing formulas, as was done first 
in the unrefined [20] and refined [lOj conifold case, and later generalized to arbitrary 
geometry without compact four-cycles in [T8| [19]. On the other hand, we will focus on 
a simpler physical derivation of BPS generating functions which uses the lift of the D- 
brane system to M-theory [23j. This also makes contact with M-theory interpretation 
of topological string theory, and allows to express BPS counting functions in terms of 
topological string amplitudes. Moreover this M-theory derivation can be extended to the 
counting of open BPS states, i.e. open D2-branes attached to additional D4-brane, which 
we are also interested in [2^ [27] . 

We start this section by reviewing the M-theory derivation of (unrefined) closed and 
open BPS generating functions. Then, to get acquainted with a crystal interpretation of 
these generating functions, we discuss their crystal interpretation in simple cases of and 
conifold. Later, using fermionic interpretation, we will generalize this crystal representation 
to a large class of toric manifolds without compact four-cycles. 

2.1 M-theory derivation 

We start by considering a system of D2 and DO-branes bound to a single D6-brane 
in type IIA string theory. It can be reinterpreted in M-theory as follows [23]. When 
additional is introduced as the eleventh dimension transversely to the D6-brane, then 
this D6-brane transforms into a geometric background of a Taub-NUT space with unit 
charge |10]. The Taub-NUT space is a circle fibration over R'^, with a circle S^j^ attaining 
a fixed radius R at infinity, and shrinking to a point in the location of the original D6-brane. 
From M-theory perspective bound states involving D2 and DO-branes are interpreted as 
M2-branes with momentum on a circle. Therefore the counting of original bound states to 
the D6-brane is reinterpreted as the counting of BPS states of M2-branes in the Taub-NUT 
space. While in general this is still a nontrivial problem, for the purpose of counting BPS 



8 



degeneracies we can take advantage of their invariance under continuous deformations of 
the Taub-NUT space, in particular under deformations of the radius R. We can therefore 
consider taking this radius to infinity, whereupon BPS counting is reinterpreted in terms 
of a gas of particles in R^. To make the problem fully tractable we have to ensure that the 
particles are non-interacting, which would be the case if moduli of the Calabi-Yau would 
be tuned so that M2-branes wrapped in various ways would have aligned central charges. 
This can be achieved when Kahler parameters of the Calabi-Yau space are tuned to zero. 
However, to avoid generation of massless states, at the same time one has to include non- 
trivial fluxes of the M-theory three-form field through the two-cycles of the Calabi-Yau 
and S)p^. In type IIA this results in the i?-field flux B through two-cycles of Calabi-Yau. 
Finally, to avoid creation of the string states arising from M5-branes wrapping four-cycles 
in Calabi-Yau, we simply restrict considerations to manifolds without compact four-cycles. 
For a state arising from D2-brane wrapping a class j3 the central charge then reads 

Z{l,^) = ^{l + B-^), (3) 

where I counts the DO-brane charge, which is taken positive to preserve the same super- 
symmetry. 

Under the above conditions, the counting of D6-D2-D0 bound states is reinterpreted in 
terms of a gas of particles arising from M2-branes wrapped on cycles ^. The excitations of 
these particles in R^, parametrized by two complex variables Zi,Z2.i are accounted for by 
the modes of the holomorphic field 

^(^i,^2) = X!"'i.'24'4'- (4) 

Decomposing the isometry group of as 5*0(4) = SU{2) x 317(2)' there are N^'"^ five- 
dimensional BPS states of intrinsic spin {m,m'). We are interested in their net number 
arising from tracing over SU{2)' spins 

m' 

The total angular momentum of a given state contributing to the index is I — li-\-l2-\-m. 
Finally, in a chamber specified by the moduli R and S, the invariant degeneracies can be 
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expressed as the trace over the corresponding Fock space 

ZbPS = {^"^ Fock<l'^° \chamber = 

= n n a - Q's^^'^-'QTncHamber 
f},m li+l2=l 
oo 

/3,m 1=1 

where the subscript chamber denotes restriction to those factors in the above product, 
which represent states which are mutually BPS 

Z{l,/3)>0 ^ g^™Q^<l. (6) 

As usual, Q = and Qs = e~^= above encode respectively the Kahler class T and the 
string coupling (we wish to distinguish carefully Qg which encodes string coupling, from 
a counting parameter q which will arise in what follows in crystal interpretation). The 
above condition on central charges is crucial in determining a particular form of the BPS 
generating functions. If we would restrict products in the formula ([5]) to factors with 
only positive /3, we would get (up to possibly some factor of MacMahon function) the 
Gopakumar-Vafa representation of the topological string amplitude. With all negative 
and positive values of /3 we would get modulus square of the topological string partition 
function. Therefore the upshot of [23] is that in general the above BPS generating function 
can be expressed in terms of the closed topological string partition function 

ZbPS = Ztop{Q)Ziop{Q ^)\chamber^ (7) 

where chamber restriction is to be understood as picking up only those factors in 
Gopakumar-Vafa product representation of Ztop for which ([6]) is satisfied. In this context 
we will often refer to the choice of a chamber as a closed BPS chamber. The (instanton 
part of the) closed topological string partition function entering the above expression is 
given by [3, 4J 

oo 

z,M) = M{qsY'' n n (1 - Q'€^T^. 

1=1 (3>0,m 

where M{qs) = Ylii^ ~ ^i) ' is the MacMahon function and x is the Euler characteristic 
of the Calabi-Yau manifold. 

To be more precise, an identification as a topological string partition function or its 
square arises if i? > in ([3]). Because R arises just as a multiplicative factor in ([3]), degen- 
eracies depend only on its sign. Therefore another extreme case corresponds to negative R 
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and B sufficiently small, when only a single D6-brane contributes to the partition function 

Z{R<0,0< B «1) = 1. (8) 

More generally, for i? < 0, BPS generating functions often (but not always) take finite 
form. 

In what follows we denote BPS generating functions in chambers with positive R by Z, 
and in chambers with negative R hj Z (and often omit the subscript BPS). Topological 
string partition functions will be denoted by Ztop, while generating functions of melting 
crystals by ordinary Z. 

The above structure can be generalized by including in the initial D6-D2-D0 configura- 
tion additional D4-branes wrapping lagrangian cycles in the internal Calabi-Yau manifold 
and extending in two space-time dimensions fl^ [2^ f27\ . For simplicity we consider a 
system with a single D4-brane wrapping a lagrangian cycle. There are now additional 
BPS states in two remaining spacetime dimensions arising from open D2-branes ending 
on these D4-branes. Their net degeneracies A^s,/?,^ are characterized, ffistly, by the 5*0(2) 
spin s whose origin is most clearly seen from the M-theory perspective [3, E]. Secondly, 
they depend on two-cycles /3 wrapped by open M2-branes, as well as one-cycles 7 on which 
these M2-branes can endj^ 

Lifting this system to M-theory we obtain a background of the form Taub-NUT x 
Calabi-Yau x S^, with the additional D4-brane promoted to M5-brane. This M5-brane 
wraps the lagrangian submanifold L inside Calabi-Yau, the time circle 5*^, and x S^js! 
inside the Taub-NUT space. A part of this lagrangian L is a torus = S*^^ x S^, which 
will lead to some modular properties of the BPS counting functions: this modularity will be 
manifest in one chamber, where the open topological string amplitude will be completed to 
the product of 6 functions. This M5-brane also breaks the 5*0(4) spatial symmetry down to 
5*0(2) X 5*0(2)'. We denote the spins associated to both 5*0(2) factors respectively by cr and 
a', and the degeneracies of particles with such spins by N^'^ . In addition to closed Kahler 
parameters Q = e~^, let us also introduce open ones related to discs wrapped by M2-branes 
z = e~°'. The real and imaginary parts of T encode respectively the sizes of two-cycles 
P and the value of the i?-field through them. The real and imaginary parts of d encode 
respectively sizes of the discs and holonomies of the gauge fields around them. Similarly 

^In case of N D4-branes wrapping the same lagrangian cycle, these states would additionally arise in 
representations R of U{N) [S]. In case of a single brane this reduces to C/(l), and such a dependence can 
be reabsorbed into a parameter specifying a choice of 7. 
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as in the closed string case, to get non-trivial ensemble of mutually supersymmetric states, 
we set the real parts of T and d to zero, and consider non-trivial imaginary parts. 

From the M-theory perspective we are interested in counting the net degeneracies of 
M2-branes ending on this M5-brane 

a' 

In the remaining three-dimensional space, in the i? — oo limit, the M2-branes ending on 
the M5-brane are represented by a gas of free particles. These particles have excitations 
in which we identify with the Zi-plane. To each such BPS particle, similarly as in the 
closed string case discussed above and in [251 HP] , we can associate a holomorphic field 

^{z,) = Y,o^iz[. (9) 
I 

The modes of this field create states with the intrinsic spin s and the orbital momentum / 
in the plane. The derivation of the BPS degeneracies relies on the identification of this 
total momentum a + 1 in the i? — oo limit, with the Kaluza-Klein modes associated to the 
rotations along S*^^ for the finite i?, following the five-dimensional discussion in P0| 

The BPS generating functions we are after are given by a trace over the Fock space 
built by the oscillators of the second quantized field and restricted to the states 

which are mutually supersymmetric. In such a trace each oscillator from Q gives rise to 
one factor of the form (1 — ql^^~^^'^ z'^)'^^ , where the exponent ±1 corresponds to the 
bosonic or fermionic character of the top component of the BPS state, 

oo 

cr,/3,7 (=1 

where the product is over either both positive or both negative 7). The parameters g, Q 
and z specify the chamber structure: the restriction to a given chamber is implemented 
by imposing the condition on a central charge, analogous to ([6]), 

ql^'-^'^Q^z^ < 1. (11) 

This condition in fact specifies a choice of both closed and open chambers. The walls of 
marginal stability between chambers correspond to subspaces where, for some oscillator, 
the above product becomes 1, and then the contribution from such an oscillator drops out 
from the BPS generating function. 
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Similarly as in the closed string case, the above degeneracies can be related to open 
topological string amplitudes, rewritten in [5j in the form 

= ( E E E ^--/^.^ -n/2j ' 

with integer Ooguri-Vafa invariants Ai'cr,^,^!^! This formula represents in fact a series of 
quantum dilogarithms 

L(z,g.) = exp (E , -n/., ) = IK^ " ^^r^'')' (^2) 

and can be written in the product form 

°° N 

zz^'iQ^ = n n n (i - q'^'^^-^") m 

cr /3,7>0n=l 

Comparing with (ITU]) we conclude that the BPS counting functions take form of the mod- 
ulus square of the open topological string amplitude 

^BPS ~ ^top {Q^^)Ztop {Q 1^ ) I chamber- (14) 

Similarly as in the closed string case, there are also a few particularly interesting cham- 
bers to consider. For example, in the extreme chamber corresponding to ImT, Imci — )■ 0, 
the trace is performed over the full Fock space and yields the modulus square of the open 
topological string partition function. In this case the quantum dilogarithms arise in pairs, 
which (using the Jacobi triple product identity) combine to the modular function Q'ijri^ in 
consequence the total BPS generating function is modular and expressed as a product of 
such functions. 

2.2 Crystal interpretation 

Closed BPS generating functions (jS]) turn out to be generating functions of statistical 
models of crystals, when parameters relevant for both interpretations are appropriately 
matched. Physical reasons for such relations have been given in [HI EH I22], and mathe- 
matical interpretation arose from works [9l [TH [15]. Such crystal interpretation arises also 

^In case of N D4-branes wrapping a lagrangian cycle this structure is again more complicated, because 
the states in R'^ arise in representations of U{N) [5]. This requires replacing the factor z""^ by the sum 
J2r Trj^y" of traces in all possible representations R of this U{N) of the matrix V encoding holonomies 
of the gauge fields. For simplicity we restrict here to the simplest case. 
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from the fermionic formulation [TB| [T7]. as we will review below. These crystals, in a more 
intricate way [27], encode also open BPS generating functions (ITU]) . However, before dis- 
cussing details of all these constructions, in this introductory section we present crystal 
models for two simplest toric Calabi-Yau manifolds, i.e. C'^ and conifold. 

is the simplest Calabi-Yau manifold. It has no compact two-cycles, so relevant BPS 
states are bound states of arbitrary number of DO-branes with a single D6-brane wrapping 
entire C'^. Their generating function is therefore expressed in terms of a single parameter 
qs = e~^\ There is just a single non-zero Gopakumar-Vafa invariant N^^q = — 1, and as 
follows from ^ this generating function coincides with the so called MacMahon function 

oo ^ 

^^^s = X{T^—^,=M{q,). (15) 

1=1 ^ 

On the other hand, the MacMahon function is a generating function of plane partitions, 
i.e. three-dimensional generalization of Young diagrams. These plane partitions represent 
the simplest three-dimensional crystal model, namely they can be identified with stacks of 
unit cubes filling the positive octant of space, as shown in fig. [TJ A unit cube located in 
position (/, J, K) can evaporate from this crystal only if all other cubes with coordinates 
{i I , j J, k < K) are already missing. A plane partition vr is weighted by the number 
of boxes it consists of |7r|, with a weight q associated to a single box, so indeed 

oo 

Z = Y^ gl-l = = 1 + g + 3g2 + 6g=^ + 13g^ + . . . = M(g), 

TT 1 = 

where p{l) is the number of plane partitions which consist of / cubes. Therefore plane 
partition generating function coincides with the BPS counting function Z = Zbps when a 
simple identification 

(ls = q (16) 

is made. From ([7]) it follows that the topological string partition function for C'^ is given 
by the square root of the MacMahon function 

Ztop = M{q,Y'\ (17) 

which is indeed true. The relevance of the MacMahon function for geometry was noticed 
for the first time in j3], and a statistical model interpretation of this result was proposed in 
[7]. We also discuss its appearance from the matrix model viewpoint at the end of section 

m 



14 



Figure 1: Plane partitions represent melting crystal configurations ofC^. 



The conifold provides another simple, yet non-trivial example of toric Calabi-Yau man- 
ifold. It consists of two C'^ patches glued into 0{—l) © 0{—l) — ?■ P^, and it has one 
Kahler class representing P^, parametrized hj Q = e^'^ . This class can be wrapped by 
D2-branes, which bind with DO-branes to an underlying D6-brane and give rise to BPS 
states in low energy theory. In this case there is already a non-trivial structure of cham- 
bers and walls, which was analyzed in |20| [2T| [HI [18]. This structure is consistent with 
M-theory derivation discussed in section [2m The generating functions of D6-D2-D0 bound 
states are parametrized by Q and qg, and therefore corresponding crystal models consist of 
two-colored three-dimensional partitions. The Kahler moduli space consists of several in- 
finite countable sets of chambers, and in each chamber relevant crystal configurations take 
form of so called pyramid partitions. These partitions are infinite or finite (respectively 
for positive and negative i? in ([2])) and their size depends on the value of the 5-field. This 
size changes discretely and the pyramid is enlarged when the value of the i?-field crosses 
integer numbers, which changes the chamber in the moduli space, as explained in previous 
subsection. Examples of such infinite pyramid partitions are given in fig. [21 and finite ones 
in fig. [21 

To write down explicitly BPS generating functions for the conifold in various chambers 
we can take advantage of their relation to the topological string amplitude dZj). The 
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Figure 2: Infinite pyramids with one and four balls in the top row, with gener- 

7Pl 

'3 



atmg junctions given respectively by Zq and Zg 



topological string partition function in this case reads 

Zt7^"''iQ) = M(g.) - QQsf, (18) 

k>l 

with the MacMahon function defined in (I15p . From this topological string partition func- 
tion we can read off Gopakumar-Vafa invariants |3l H] 

N"^=o = -2, ATj^ii = 1. 

Using the relation ([7]) we can now present conifold closed BPS generating functions 
in several sets of chambers. In the first set of chambers we consider R > and positive 
B &]n,n+ 1[ (for n > 0). Firstly, for small B, there is so-called non-commutative chamber 
discussed first by Szendroi [14J, which corresponds to ri = 0. In this case the pyramid 
crystal has just a single ball in the top row, as in the left panel in fig. [21 and the BPS 
generating function is given by the square of the topological amplitude. On the other 
hand, for large B, i.e. n — t- oo, we reach commutative chamber in which the length of the 
top row extends to infinity. In this case the BPS generating function agrees, up to a single 
factor of MacMahon function, with the topological string amplitude. In between there are 
chambers with n -\- 1 balls in the top row, for which 

k>l k>n+l 

These BPS generating functions are related to pyramid generating functions with two 
colors go and gi upon the identification (which generalizes (fT6|l in case) 

^conifold chambers : = goQi, Q = -g^gi- 
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Indeed, with this identification the above counting functions agree with those of two-colored 
pyramid crystals with n + 1 yellow balls in its top row 

^r™"^^'(go, qi) = M{q,q,f J] (1 + go'g^')'"" + qo^t^f^^- (20) 

fc>n+l k>l 

In the second set of chambers we have < and positive B E\n — l,n[ (for n>l). It 
extends between the core region with a single D6-brane ([S]) and the chamber characterized 
by so-called Pandharipande-Thomas invariants (for the flopped geometry, or equivalently 
for anti-M2-branes) . The BPS generating functions read 

n— 1 j 

^conifold = JJ (1 - '^y. (21) 

The corresponding statistical models were shown in |16[ EH [IB] to correspond to finite 
pyramids with n — 1 stones in the top row, as shown in figure [3l In this case the generating 
functions of such partitions are equal to 

n-l 

Zvyr'^rma^q,^ ^i) = + qt' qr'^ ■ (22) 

The equality ^con.j/oid ^ ^pyramid g^j^jggg upon an identification 

^comfoid chambers : q-^ = qoqi, Q = -q^qi. 

There are two other sets of chambers characterized by the negative value of the S-field, 
for which BPS generating functions are completely analogous to those given above. 

Above we presented just the simplest examples of crystal models. Using fermionic 
formulation presented below one can find other crystal models for arbitrary toric geometry 
without compact four-cycles. Let us also mention that those models can be equivalently 
expressed in terms dimers. In particular the operation of enlarging the crystal, as in the 
conifold pyramids, corresponds to so called dimer shuffling [15j. Dimers are also closely 
related to a formulation using quivers and associated potentials, which underlies physical 
derivations in [2T| [22]. 

3 A little background — free fermions and matrix models 

In this section we introduce some mathematical background on which the main results 
presented in this review rely. In section 13.11 we start with a brief presentation of toric 
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Figure 3: Finite pyramids with m = 1,2,3 stones in the top row (respectively 
left, middle and right), whose generating functions are given by Z^^l^^'^ (note 
that z'^'^'^'^^'^ = 1 corresponds to an empty pyramid corresponding to the pure 
D6-hrane). 



Calabi-Yau manifolds and introduce the notation which we use in what follows. In section 
13.21 we introduce free fermion formalism. In section 13.31 we introduce basics of matrix 
model formalism. Our presentation is necessarily brief and for more detailed introduction 
we recommend many excellent reviews on each of those topics. 

3.1 Toric Calabi-Yau three- folds 

Some introductory material on toric Calabi-Yau manifolds, from the perspective rel- 
evant for mirror symmetry and topological string thoery, can be found e.g. in [39j. In 
this section our presentation is brief and mainly sets up the notation. Toric Calabi-Yau 
three- folds arise as the quotient of C^'^, possibly with a discrete set of points deleted, by 
the action of (C*)'' with certain weights. The simplest toric three-fold is C^, which corre- 
sponds to the trivial choice k = 0. The resolved conifold, which we already discussed in 
section [2l2| corresponds to k = 1 and a choice of weights (1, 1, —1, —1), which represent a 
local bundle 0{—l)®0{—l) — J-P^. The structure of each toric three- fold can be encoded 
in a two-dimensional diagram built from trivalent vertices. Finite intervals joining two 
adjacent vertices represent local neighborhood inside the manifold. Equivalently one 
can consider dual graphs. Examples of toric diagrams and their duals for C'^, conifold and 
resolution of C'^/Z2 singularity are given in fig. H] (the notation F-i- at each vertex will be 
explained in what follows). 

A closed loop in a toric diagram represents a compact four-cycle in the geometry. As 
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Figure 4: Toric graphs for C^, conifold and resolution of C^/Z2. 



follows from the reasoning in section 12. ![ in the context of BPS counting we are forced to 
restrict considerations to manifolds which do not have such four-cycles. Apart from a few 
special cases, there is an infinite class of such geometries whose dual diagrams arise from 
a triangulation, into triangles of area 1/2, of a long rectangle or a strip of height 1. A 
toric diagram arises as a dual graph to such a triangulation. From each vertex in such a 
toric diagram one vertical line extends to infinity and crosses either the upper or the lower 
edge of the strip. Two such consecutive lines can emanate either in the same or in the 
opposite direction, respectively when they are the endpoints of an interval representing 
with local 0{—2) © (9 or 0{—l) © 0{—l) neighborhood. Example of a generic diagram of 
this kind is shown in fig. [8] below. 

Let us denote independent P^'s, starting from the left end of the strip, from 1 to A^, 
and introduce corresponding Kahler parameters Qi = e~'^\ i = 1, . . . , N. Moreover to each 
toric vertex we associate a type ti = ±1, so that tj+i = ti if the local neighborhood of P^ 
(represented by an interval between vertices i and i + is 0{—2)(B0; if this neighborhood 
is of 0{—l) © 0{—l) type, then tj+i = — tj. The type of the first vertex we fix as ti = +1. 
In figures m and [H] these types are denoted by © and ©. The types tj will be used much in 
the construction of fermionic states in section 14.21 

As explained in section [2rT| the BPS generating functions can be expressed in terms of 
(the instanton part) of topological string amplitudes. For the above class of geometries, 
arising from a triangulation of a strip, these amplitudes read 

ZtoMd = M{q,)^ H H (l-qi {Q,Q,+^ ■ ■ ■ Q,_i)) ' ' . (23) 

1=1 l<i<j<N+l 
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3.2 Free fermion formalism 



Formalism of free fermions in two dimensions is well known |^ HI] and ubiquitous in 
literature on topological strings and crystal melting jTJ [151 EZl US SH]. The main purpose 
of this section is therefore to set up the notation which we will follow in the remaining 
parts of this paper. 

The states in the free fermion Fock space are created by the (anti-commuting) modes 
of the fermion field 

^(^) = ^i^k+\/2Z~^~^ : 1p*{z) = ^1p*k+l/2Z'^'^, {^fc+l/2,Ci-l/2} = ^k,l 

kez fcgz 
on the vacuum state |0). There is one-to-one map between such fermionic states 

d 

1/^) = n^-«-i/2^-^-i/2|0)' with ai = fii-i, bi = fil-i, 

1=1 

and two-dimensional partitions fi = (yUi,/i2, • • •A*;)' shown in fig. 13 The modes am of 
the bosonized field dcj) =: ■ip{z)ip*{z) : satisfy the Heisenberg algebra [am,tt-n] = n6m,n- 



I 
I 




Figure 5: Relation between Young diagrams and states in the Fermi sea. 



We introduce vertex operators 

T^{x) = e^">o "^"i", T'^{x) = e^">o « °±", (24) 

which act on fermionic states \fi) corresponding to partitions fi as |13| HH [15] 

r-(^)l/i) = $^a:I^H^I|A), T4xM = a;!-!"!^! |A), (25) 

r'_(x)|/.) = Yl ^'""'"'|A), r'^ixM = J2 ^'^'-'"|A). (26) 
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The interlacing relation -< between partitions is defined as 



^ <^ Ai > /ii > A2 > > A3 > . . . . (27) 

The operator F' is the inverse of F with negative argument. These operators satisfy 
commutation relations 

F+(x)F_(i/) = ^—T4y)T4x), (28) 
I — xy 



T'4x)T'_iy) = ^—-r_iy)T'4x), (29) 



1 

1 — xy 

T'4x)T4y) = {l + xy)T4y)T'4x), (30) 
F+(x)F'_(y) = il + xy)T'_{y)T+{x). (31) 

We also introduce various colors qg and the corresponding operators Qg (a hat is to 
distinguish them from Kahler parameters Qi) 

Q.|A) = gl^l|A). (32) 

These operators commute with vertex operators up to rescaling of their arguments 

^+{x)Qg = QgT+{xqg), ^'+{x)Qg = QgT'^ixqg), (33) 

QgT.ix) = r_{xqg)Qg, Qg^'-{x) = r'_{xqg)Qg. (34) 



3.3 Matrix models 

In matrix model theory, or theory of random matrices, one is interested in properties 
of various ensembles of matrices. Excellent reviews of random matrix theory can be found 
e.g. in |46] or, in particular in the context of topological string theory, in |47| . In matrix 
model theory one typically considers partition functions of the form 

Z= f VUlle~i^^'''^''\ (35) 

where V = V{U) isa matrix potential, and VU is a measure over a set of matrices of interest 
U of size A^. Typically it is not possible to perform the above integral, however special 
techniques allow to determine its formal 1/A^ expansion. These techniques culminated 
with the formalism of the topological expansion of Eynard and Orantin [38] which, in 
principle, allows to determine entire 1/A^ expansion of the partition function recursively. 
This solution is determined by the behavior of matrix eigenvalues, whose distribution 
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among the minima of the potential, in the continuum hmit, determines one- dimensional 
complex curve, so-called spectral curve. The spectral curve is also encoded in the leading 
expansion of the so-called resolvent, which is defined as the expectation value u{x) = 
(Tr^^) computed with respect to the measure (l35l) . 

In the context of BPS counting and topological strings, unitary ensembles of matrices 
of infinite size arise. In this case the matrix model simplifies to the integral over eigenvalues 
Ua, with a measure which takes form of the unitary Vandermonde determinant 



VU = l[duo. Yl 



2 

5 



The issue of infinite matrices is a little subtle, however it can be taken care of by considering 
matrices of large but finite size N, and subsequently taking N ^ oo limit. For finite one 
can find the resolvent, and in consequence the spectral curve, using a standard technique 
of so-called Migdal integral. This requires redefining V to the standard Vandermonde form 
[m [in], as well as introducing 't Hooft coupling T 

V + T\ogz, T = Ng,. (36) 

The form of the Migdal integral depends on the number of cuts into which eigenvalues 
condense in large N limit, and this number of cuts determines the genus of the spectral 
curve. In our context only single-cut situations will arise, for which the spectral curve has 
genus zero. In this case the Migdal integral determines the resolvent as 

I fdz a.Vjz) V(p - a)(p - b) 

SO that the integration contour encircles counter-clockwise the endpoints of the cut a and 
b. A proper asymptotic behavior of the resolvent is imposed by the condition 

lim ijj{p) = -. (38) 

p^oo p 

Then the spectral curve is determined as a surface on which the resolvent is unambiguously 
defined, i.e. it is given by an (exponential) rational equation automatically satisfied by 
p and u){p). There is also an important consistency condition for the resolvent: when 
computed on the opposite sides of the cut uj{p)±, it is related to the potential as 

a;+(p)+c._(p) = ^^. (39) 
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On the other hand, a difference of these values of the resolvent on both sides of the cut 
provides eigenvalue density 

p{p) =uj+{p) -uj-{p). (40) 

It has been observed in several contexts that topological strings on toric manifolds can 
be related to matrix models, whose spectral curves take form of the so-called mirror curves. 
Mirror curves arise for manifolds which are mirror to toric Calabi-Yau manifolds |39[ 
mj. For toric manifolds, their mirror manifolds are determined by the following equation 
embedded in four-dimensional complex space 

Z1Z2 = H{x,y). 

The mirror curve is the zero locus of H{x, y), i.e. it is given as H{x, y) = 0. More precisely, 
X, y are C* variables, and it is often convenient to represent them in the exponential form 
X = u", y = e"", with u,v ^ C For example, for C'^ and the conifold they take the following 
form 

Hc3{x,y) = X + y + xy = 0, Hconifoid{x,y) = x + y + xy + Qx^ = 0, (41) 

where Q encodes the Kahler parameter of the conifold. Schematically mirror curves arise 
from thickening edges of the toric graphs, as shown in fig. |6l 

One of the first relations between topological strings for toric manifolds and matrix 
models were encountered in [ ^ l50 ] . where it was shown that the spectral curve of a unitary 
matrix model with a gaussian (i.e. quadratic) potential agrees with the above mirror curve 
Hconifoid{x,y) = in (HT]) . with 't Hooft coupling T = QgN encoded in Q = e^^ . At 
the same time it was shown that the matrix model partition function reproduces the 
topological string partition function. More recently these ideas became important in view 
of the remodeling conjecture [5T| [57] . which states that the solution to loop equations 
in the form found by Eynard and Orantin [l8j, applied to the mirror curve, reproduces 
topological string partition functions. The method of [48j works for arbitrary curves, not 
necessarily originating from matrix models. Nonetheless, it is indeed possible to construct 
matrix models whose partition functions do reproduce topological string amplitudes, and 
whose spectral curves coincide with appropriate mirror curves [29| [30| |52| [53| |55| [56] . 

One of our aims is to provide matrix model interpretation of BPS counting. It is natural 
to expect such an interpretation in view of an intimate relation between BPS counting 
and topological string theory discussed in section 12. H and the above mentioned relations 
between topological strings and matrix models. As we will see in what follows, there are 
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indeed unitary matrix models which naturally arise in the context of BPS counting and 
its fermionic formulation. Among the others, our task will be to analyze them using the 
above mentioned Migdal method. 




Figure 6: Toric diagrams for and conifold and corresponding mirror curves. 



4 Fermionic formulation of BPS counting functions 

Having introduced all the ingredients above, we are now ready to present fermionic 
formulation of BPS counting. To start with, in section 14.11 we present the idea of such a 
formulation in the simplest example of C'^. In section we introduce a general fermionic 
formalism, and in section 14.31 we provide its crystal interpretation. We illustrate the use 
of our formalism in section H741 revisiting C'^ example, as well as in explicit case of C'^/Z^r, 
and conifold geometry. 

4.1 The idea and example 

As explained in section \2.2\ the generating function of bound states of DO-branes to a 
single D6-brane is given by the MacMahon function, and the corresponding crystal model 
takes form of the counting of plane partitions [7j. Let us slice each such plane partition 
by a set of parallel planes, as shown in fig. [71 In this way on each slice we obtain a 
two-dimensional partition /i, and it is not hard to see that each two neighboring partitions 
satisfy the interlacing condition (127]) . Recalling that such a condition arises if we apply 
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r±(l) operators f l25p to partition states, we conclude that a set of all plane partitions can 
be built, slice by slice, by acting with infinite sequence of r-|-(l) on the vacuum. To count 
each slice /i with appropriate weight g''^' we also need to apply weight operator Q defined 
in (132|) . Therefore the generating function of plane partitions can be represented as follows 

z = = (o|...Qr+(i)Qr+(i)Qr+(i) I Qr_(i)gr_(i)Qr_(i)g...|o) = 

= (0| . . .r+(g2)r4g)r4l)r_(g)r_(g2)r_(g=^) . . . |0) = (42) 

oo ^ 

= n 1 _ ,....-1 = ^(^)- 

h,h=l ^ 

In the first line we implicitly introduced two states and |^^-), defined by an infinite 
sequence of r_|_ (respectively r_) operators, interlaced with weight operators Q and acting 
on the vacuum. To confirm that this correlator indeed reproduces the MacMahon function, 
the second line can be reduced to the final infinite product using commutation relations 
(128|) and f lM|) . We can also represent insertions of r±(l) operators graphically by arrows, 
so that the above computation can represented as in fig. [7] (right). 




Figure 7: Slicing of a plane partition (left) into a sequence of interlacing two- 
dimensional partitions (right). A sequence ofr± operators in which create 
two-dimensional partitions is represented by arrows inserted along two axes. 
Directions of arrows — )■ represent interlacing condition >- on partitions. We 
reconsider this example from a new viewpoint in figure ITU. 
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In what follows we present a formalism which allows to generalize this computation to 
a large class of chambers, for arbitrary toric geometry without compact four-cycles. 

4.2 Toric geometry and quantization 

We wish to reformulate BPS counting in the fermionic language in a way in which we 
associate to each toric manifold a fermionic state, such that the BPS generating function 
can be expressed as an overlap of two such states, generalizing case ( H2|) . At the 
same time the construction of such a fermionic state is supposed to encode the structure 
of the underlying crystal model (generalizing plane partitions in fig. [7]). An important 
difference between and other geometries is the existence of many Kahler moduli and 
correspondingly many chambers, for which BPS generating functions change according to 
wall-crossing formulas. To take care of these changes in the fermionic formalism we need 
to introduce special wall- crossing operators. 

4.2.1 Toric geometry and fermionic operators 

In what follows we use the notation introduced in section 13.11 in particular to each 
vertex of the toric diagram we associate its type ti = ±1, see also fig. [HI We start with a 
construction of fermionic states associated to a given toric Calabi-Yau manifold (without 
compact four-cycles). First we need to introduce several operators which are building 
blocks of such states. The structure of these operators is encoded in the toric diagram of 
a given manifold. Namely, these operators are given by a string of + 1 vertex operators 
r^(x) (defined in (IMl) ) which are associated to the vertices of the toric diagram; the type 
ti determines the type of a vertex operator as 

Tl=+\x) = r±(x), rfix) = r'i(x). (43) 

In addition the string of operators r^(a;) is interlaced with A^ + 1 operators Qi representing 
colors qi, for i = 0,1, . . . , N. Operators Qi, . . . ,Qn are associated to in the toric 
diagram, and there is an additional Qq. We also define 

Q = QoQi ■ ■ ■ Qn, q = qoqi ■■■qN- (44) 

Therefore the upper indices of T%{x) and a choice of colors of the operators which we 
introduce below are specified by the data of a given toric manifold. As we will see, a 
sequence of lower indices ± is determined by the chamber we are going to consider. 
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Figure 8: Toric Calabi- Yau manifolds represented by a triangulation of a strip. 
There are N independent 's with Kdhler parameters Qi = e~'^\ and + 1 
vertices to which we associate T and T' operators represented respectively by © 
and signs. Yellow intervals, which connect vertices with opposite signs, rep- 
resent 0{—l) © 0{—l) — i- local neighborhoods. Red intervals, which connect 
vertices with the same signs, represent 0{—2) © (9 — t- P^ local neighborhoods. 
The first vertex on the left is chosen to be ®. 



Now we can associate several operators to a given toric manifold. Firstly we define 

A±{x) = Tl{x)Q^Tt{x)Q2 ■ ■ ■ Tt{x)QNT'r\x)Qo. (45) 
Commuting all QiS using fl34p we also define the following operators 

A^{x) = A^{x) = F^j (xg)r^ {^)T^^ (^) ■ ■ ■ F^^ ( ), (46) 

qi qiq2 qiq2---qN 

A_{x) = A_{x) = F*_^(x)F*_2(xgi)F*_^(xgig2) ■ ■ ■ F*_^+^ (xgigag^v)- (47) 
In addition, we define the above mentioned wall-crossing operators 

Wp{x) = (^r'l{x)Q,r'^{x)Q2 - ■ ■T'l{x)Q,^(^r'^+\x)Q,+, - ■ ^ (48) 
W'pix) = (^T%ix)QiT'_^ix)Q2---T'^{x)Q,)(^T'l+\x)Qp^^ (49) 

Here the order of F and F' is the same as for A± operators, and the difference is that now 
there are subscripts =]= on first p operators and ± on the remaining ones. 
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We often use a simplified notation when the argument of the above operators is x = 1 

Z± = A±(i), /i± = A±(i), Wp = Wp(i), w'^ = w'^{i). 

4.2.2 Fermionic formulation and quantization 

Above we associated operators A± to each toric geometry with a strip-hke toric diagram. 
From these operators we can build the following states in the Hilbert space of a free fermion 

n 

|fi±) G (50) 

which we define as follows 

(fi+l = (0|...A+(l)Z+(l)A+(l) = (0|...A+(g2)A+(g)A+(l), (51) 
= Z_(l)A_(l)A_(l)...|0)=A_(l)A_(g)A.(g2)...|0). (52) 

These states encode the full instanton part of the topological string amplitudes. Namely, 
as shown in [16j, 

z = {n+\n_) (53) 

is equal to the BPS partition function Z in the non- commutative chamber 

Z = Z = \Ztap\^ = Ztop{Qi)Ztop{Qi (54) 

where Ztop{Qi) is given in (123|) . The above equality holds under the following identification 
between qi parameters (which enter the definition of |fi±)) and physical parameters Qi = 

e~^^ and = e~^°: 

Qi = {titi+i)Qi, qs = q = Qoqi ■■■qN- (55) 

We will provide a proof of ( 153|) in section iG.l.ll in a more general setting of refined invariants. 

The states have non-trivial structure and encode the information about the non- 
commutative chamber. It turns out that the fermionic vacuum |0) itself also encodes some 
interesting information. We recall that there is another extreme chamber representing just 
a single BPS state represented by the D6-brane with no other branes bound to it. This 
multiplicity 1 can be understood as 

Z = Z={0\0) = 1, (56) 

and as we will see below, starting from this expression we can use wall-crossing operators 
to construct BPS generating functions in an infinite family of other chambers. 
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4.2.3 Other chambers and wall-crossing operators 

In the previous section we associated to toric manifolds the states \Q±), whose overlap 
reproduces the BPS generating function in the non-commutative chamber (153|) . Now we 
wish to extend this formalism to other chambers. As discussed in section 12. ![ in a given 
chamber, the allowed bound states we wish to count must have positive central charge 

Z{R,B) = hn + (3-B)>0. 
K 

Firstly, the information about R and B must be encoded in the fermionic states which we 
wish to construct. It turns out that the choice of positive or negative R is encoded in the 
choice of the ground state 

R>0 — ¥ |fi±), R<0 — ¥ |0), (57) 

which generalizes the extreme cases (155]) and f l56p . 

On the other hand, the value of the field B is encoded in the insertion of additional wall- 
crossing operators, such as those defined in (H8|) and ( H9|) . In particular these two types 
of operators are sufficient if we wish to consider only these chambers, which correspond 
to a flux of the 5-field through only one, but arbitrary in the manifold. For simplicity 
below we consider only this set of chambers. Denoting this as p, it can be shown that 
insertion of n copies of operators Wp or W'^ creates respectively n positive or negative 
quanta of the flux through p'th P^. 

Therefore, schematically, the generating functions in chambers with R > read 

Zn = {n+\{Wr\n^), (58) 

and those with R < read 

^n=(0|(W)"|0), (59) 

with appropriate form of wall-crossing operators. More precisely, depending on the signs 
of R and B, we need to consider four possible situations, which we present below. The 
proofs of all statements below, corresponding to these four situations, can be found in [,16j. 

• Chambers with R <0, B > 

Consider a chamber characterized by positive R and positive i?-field through p'th 
two-cycle 

R<0, Be]n-l,n[ for 1 < n G Z. 
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The BPS partition function in this chamber contains only those factors which include 
Qp and it reads 



n-l p N+l 



^».=nn n o 



1. 



■ 1 1 ,1 QsQs+l ■ ■ ■ Qr-l 

1=1 s=l r=p+l 

This can be expressed as the expectation value of n wall-crossing operators Wp 

^n|p=(0|(Wp)"|0)=Z„,p, (60) 

under the following identification of variables 

Qp = {tptp+i)qpq'^, Qi = {titi+i)qi for i ^ ^ ^' ^^"^^ 

A special case of this result is the trivial generating function (156 p representing a 
single D6-brane. 

Chambers with i? > 0, 5 > 

In the second case we consider the positive value of R and the positive flux through 
p'th 

i?>0, Be]n,n + l[ for < n G Z. 

Denote the BPS partition function in this chamber by Zn\p. We find that the ex- 
pectation value of n wall-crossing operators Wp in the background of \Vt) has the 
form 

Z.I, = {n^\{Wpr\n^) = M(l,g)^+^zi'Jzi;)z(^, (62) 

where Z^^^ does not contain any factors (g^ ■ ■ ■ gr-i)^^ which would include g,, while 
Z^jp contains all factors qs - ■ ■ qr-i which do include g,, and Z^^^ contains all factors 
(q's ■ ■ ■ Qr-i)'^ which also include g,: 

= fi n (i-(M.)--^^— )"''''^Yi-(Ms)gw+i---g.-i'"*^*^' 

, -, , ^ qsqs+i qr—i' ^ ' 

'=1 p^s,r+l(Zl,N+l 



^n\p = n n (i-(M.)g'+"g.g.+i---gr-i 

'=1 pes,r+lCl,N+l 

= n n i-(w-) — - — 

l=n+l p£s,r+lCl,N+l 
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We see that the identification of variables 

Qp = (tptp+i)qpq'^, Qi = (titi+i)qi for i p, Qs = q (63) 

reproduces the BPS partition function 

Zn\p = Zn\p. (64) 

When no wall-crossing operator is inserted the change of variables reduces to fl55|) and 
we get the non-commutative Donaldson-Thomas partition function (15^ . Zq\p = Z. 

• Chambers with R < 0, B < 

Now we consider negative R and negative 5-field 

R<0, B e]-n-l,-n[ for < n G Z. 

For such a chamber the BPS partition function reads 

i=l s=l r=p+l 

Now we find the the expectation value of n wall-crossing operators is equal to 

z'^\p = mw'^r\Q) = z'^\p, (65) 

under the change of variables 

Qp = {tptp+i)qpqj"', Qi = {titi+i)qi for i ^ p, ^ ^' ^^^^ 

Now an insertion of Wp has an interpretation of turning on a negative quantum of 
S-field, and the redefinition of Qp can be interpreted as effectively reducing tp by 
one unit of Qg. As already discussed 

Z'o\p = (0|0) = 1 

represents a chamber with a single D6-brane and no other branes bound to it. 

• Chambers with R > 0, B < 

In the last case we consider positive R and negative B 

R>0, > B e] -n,-n + l[ for 1 < n G Z. 
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We denote the BPS partition function in this chamber by Z',^^^. We find that the 
expectation value of n operators w'^ in the background of \^±) has the form 

z\ = = M(i, qf^' z'^^ (6?) 

where Z'^^_^ does not contain any factors (g^ ■ ■ -Qr-i)^^ which would include qp, Z'^^ 
contains all factors gs ■ ■ ■ qr-i which do include qp, and Z^^,J contains all factors 
■ ■ ■ Qr-i)^^ which also include qp-. 



= n n — )"^7i-(M.)gw+i---9.-i 

, , , ^ ysys+l ■ ■ ■ yr— 1 ' ^ 

'=1 p^s,r+lCl,N+l 
oo 



-1 



'pes,r+lCl,Af+l 



7'(2) _ TT TT h _ {+ + \ 1 



c = n n 



, , qsQs+l' ' ' Qr-l 

'=1 pes,r+lCl,JV+l 

Under the change of variables 

Qp = {tptp+i)qpq~'^~^ , Qi = {titi+i)qi for i p, Qs = Q- (68) 

this reproduces the BPS partition function 

K\p = Z'nip- (69) 

We note that both Z'^^^ with the above change of variables, as well as Zq\p given in (1^2]) 
with a different change of variables in f l63|) . lead to the same BPS generating function 
Z which corresponds to the non-commutative Donaldson-Thomas invariants. 



4.3 Crystal melting interpretation 

In the previous section we found a free fermion representation of D6-D2-D0 generating 
functions. The fermionic correlators which reproduce BPS generating functions automat- 
ically provide melting crystal interpretation of these functions [16j, generalizing models of 
plane partitions (for C^) or pyramid partitions (for the conifold), presented in section [2l2l 
These crystals are also equivalent to those found in |17[ 12^ . 

The crystal interpretation is a consequence of the fact that all operators used in the 
construction of states \^l±), as well as the wall-crossing operators, are built just from vertex 
operators F-i- and F-i- with argument 1, and color operators Qj. As follows from (125|) and 



32 



(I26p . insertion of these vertex operators is equivalent to the insertion of two-dimensional 
partitions satisfying interlacing, or transposed interlacing conditions. An (infinite) se- 
quence of such interlacing partitions effectively builds up a three-dimensional crystal. A 
relative position of two adjacent slices is determined by a type of two corresponding ver- 
tex operators. On the other hand, insertions of color operators have an interpretation 
of coloring the crystal. The colors Qi appear in the same order in each composite oper- 
ator, so these colors are always repeated periodically in the full correlators. Therefore, 
three-dimensional crystals are built of interlacing, periodically colored slices. 

To get more insight about a geometric structure of a crystal it is convenient to introduce 
the following graphical representation. We associate various arrows to the vertex operators, 
as shown in fig. [9l These arrows follow the order of the vertex operators in the fermionic 
correlators, and are drawn from left to right, or up to down (either of these directions is 
independent of the orientation of the arrow). Following the order of the vertex operators 
in a given correlator, and drawing a new arrow at the end of the previous one, produces a 
zig-zag path which represents a shape of the crystal. The coloring of the crystal is taken 
care of by keeping track of the order of Qi operators, and by drawing at the endpoint 
of each arrow a (dashed) line, rotated by 45°, colored according to Qi which we come 
across. These lines represent two-dimensional slices in appropriate colors. In this way the 
corners of two-dimensional partitions arising from slicing of the crystal are located at the 
end-points of the arrows. The orientation of arrows represents the interlacing condition 
(i.e. arrows point from a larger to smaller partition). The interlacing pattern between two 
consecutive slices corresponds to the types of two consecutive arrows. Finally, the points 
from which two arrows point outwards represent those stones in the crystal, which can be 
removed from the initial, full crystal configuration. In fermionic correlators these points 
correspond to followed by F^ operators. We illustrate this graphical construction in a 
few examples in the next section. 

4.4 Examples 
4.4.1 Revisiting 

Let us reconsider geometry which motivated our discussion in section 14.11 In this 
case, the dual toric diagram consists just of one triangle, see fig. (left), so there is just 
one vertex and only one color Qq = Q, and the operators (145|) take form 

A± = F±(i)g. 
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r+ r; r_ 



Figure 9: Assignment of arrows. 



In consequence the BPS partition function fl53p takes exactly the form fl42p . 

The crystal structure can be read off from a sequence of arrows associated to A± 
operators, following the rules in figure [9l This gives rise to the crystal shown in fig. [10] 
(right). This is the same crystal as in fig. [TJ which represents plane partitions, however 
now seen from the opposite side. 

4.4.2 Orbifolds CVZ^+i 

Now we consider the resolution of C'^/Z^^^ orbifold. In this case the toric diagram 
takes form of a triangle of area (A^ + l)/2, see fig. [11] (left). There are independent P^'s 
and + 1 vertices of the same ti = +1, and operators in f l45p take the form 

Z± = r±(i)Qir±(i)g2 • • • r±(i)Q^r±(i)go. 

In the non-commutative chamber the corresponding crystal consists of plane partitions, 
however with slices colored periodically in + 1 colors. The partition function in the 
non-commutative chamber is given by (153|) . 

If we turn on an arbitrary 5-field through a fixed P^, the structure of wall-crossing 
operators gives rise to modified containers, see e.g. fig. [TT] (middle) . In particular enlarging 
the 5-field by one unit adds one more yellow corner to the crystal. 

The crystals corresponding to i? < are also easy to find. In the extreme chamber 
we get a trivial (empty) crystal, representing a single D6-brane ( !56|) . Adding wall-crossing 
operators results in a crystal with several corners, finite along two axis (and extending 
infinitely along the third axis), as shown in fig. [11] (right). 
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Figure 10: Toric diagram for (l^ft) consists of one © vertex. Operators A±. 
involve a single r±, and have a simple arrow (lower left), as follows from figure 
0. The correlator is translated into a sequence of arrows, with rotated 
dashed lines representing insertions of interlacing two-dimensional partitions. 
The resulting figure (right) represents plane partitions crystal model, the same 
as in figurcY^ hut now seen from the bottom. 



4.4.3 Resolved conifold 

We already presented pyramid crystals for the conifold in section 12.21 They arise 
from our formalism as follows. The dual toric diagram for the conifold, see fig. [12] (left), 
consists of two triangles and encodes a single (A^ = 1) P^. Two vertices of the toric diagram 
correspond to two colors Qi and Qo? so that 

Q = QiQo, q = qiqo- 

The operators (H5|) in this case read 

A±{x) = r±(x)Qir^(x)Qo, 

while (HH]) and (H?]) are 

A+{x) = T+{xq)T'^{xq/qi), A_{x) = r_(x)r_(xgi), 

and they satisfy 

A f \ A f \ i^ + ^yq/qi)i'^ + ^yqqi) A f \ a f \ 

A+{x)A4y) = _ ^y^y A4y)A+{x). 
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Figure 11: Toric diagram for the resolution of C^/Z^'*'^ geometry has N + 
1 vertices of the same type ©. Left: toric diagram for N = 1 and arrow 
representation of and Wi. In the non- commutative chamber this leads to 
the same plane partition crystal as in fig. [73, however colored now in yellow and 
red. Middle: for the chamber with positive R and 2 < B < 3 the crystal develops 
two additional corners and its generating function reads Z2\\ = 
Right: for negative R and positive n — 1 < B < n the crystal is finite along 
two axes and develops n — 1 yellow corners; its generating function for the case 
of n = 5 shown in the picture reads Z^^^i = (0|(H^i)^|0) (two external arrows, 
corresponding to r_ acting on (0| and r_|_ acting on |0), are suppressed.) 



The quantum states (l5Ti) and (152|) take form 



A_(l)A_(g)A_(g2)...|0), 
{0\...A^iq')A4q)AM 



(70) 
(71) 
(72) 



and the wall-insertion operators f HHj) f H9|) are 



Wi{x) = T4x)QiT'4x)Qo 



Wiix) = T+{x)Q^T'_{x)Qo. 



(73) 



Therefore the fermionic correlators take form 




{n+\{w^r\n^) 
(o|(Wi)"|o). 



(74) 
(75) 
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Figure 12: Left: toric diagram for the conifold and arrow representation 
and Wi. Right: for chambers with negative R and positive n — 1 < B < n the 
crystals are given by finite pyramid partitions with n — 1 additional corners, 
represented by n — 1 stones in the top row (the figure shows the case n = 4). 
The generating function is given by Zn\i = (0|(iyi)"|0) which reproduces the 
result [2^1 . 



and encode generating functions (120|) and f l22|) introduced in section 12.21 In the non- 
commutative chamber we get the result found first in |14j, Zq\i = while a single 
D6-brane is encoded in Z = (0|0) = 1. These crystals are shown in fig. [T^ (right) and[T31 

5 Matrix models and open BPS generating functions 

In this section we explain how matrix model formalism can be applied to analyze BPS 
counting functions. In the first part, 15. 1[ we explain how to relate fermionic formalism, 
derived in the previous section, to matrix model representation. In sect ion [5l2| we illustrate 
how to construct matrix models for the closed non-commutative chamber. In subsection 
15.31 we analyze in detail BPS generating functions for the conifold for all chambers with 
-R > 0, and derive corresponding spectral curves. We discuss how these curves relate to 
(and generalize) mirror curves, which we find (as we should) in the commutative chamber. 



37 



//////////// 
//////////// 

Figure 13: Conifold crystal in the chamber with positive R and 2 < B < 3 
takes form of pyramid partitions with 3 stones in the top row. Its generating 
function is given by Z2\i = {Q+\{Wi)'^\Q-) . 



In subsection 15.41 we reveal that matrix model representation in fact encodes open BPS 
generating functions, which can be identified with matrix model integrands. 

5.1 Matrix models from free fermions 

Let us explain how to relate fermionic representation of BPS amplitudes, introduced 
in section 14. 2[ to matrix models. This relies on introducing into fermionic correlators 
representing BPS generating functions, such as ( 15 3 p or (1^2]) . a special representation of 
the identity operator I. The representation we are interested in also consists of infinite 
product of vertex operators and arises as follows |29]. Firstly, we can use the representation 
as a complete set of states I = \R){R\, which represent two-dimensional partitions. Using 
orthogonality relations of U{oo) characters xr^ ^"^^ the fact that these characters are given 



38 



in terms of Schur functions xr = Sr{z) for z = {zi, Z2, ^3, . . ■), we can write 

I = 5^|i?)(i?| = 5^W|P)(i?| = 

R P,R 



spt{z)sRtiz)\P){R\ = 

P,R 

= /"pf/(nr-(^J|o))((oinr'+fe^)). (76) 

a a 

When such a representation of the identity operator is introduced into (153|) or (162|) (or any 
other correlator of similar structure) we can commute away T'^^ operators and get rid of 
operator expressions. For example, inserting the above identity operator in the string of 
operators in (155]) . leads to a matrix model with the unitary measure 

00 fc— 1 

= (oin^+(i)iiin^+wi^i^-)= 

i=k j=0 

00 



vu{o\ n ^+(i)inr-(^a)io)(oinr;(.-^)in^+(i)i^i^-)= 

i=k+l a a j=0 

f!;{q,Q,) /l?f/ne-^^"(^"^ (77) 



a 



The product over a represents distinct eigenvalues z^- Note that we have inserted I at 
the position k in the string of v4+(l) operators. In particular this affects the form of 
the resulting potential V^{z). Moreover, apart from matrix integral, we find some overall 
factors fn{q,Qi) which take form of various infinite products. They arise, in a generic 
chamber, from commutations between F-i- ingredients of wall-crossing operators, and 
ingredients of |i7zp) states. In the closed non-commutative chamber n = these factors are 
trivial, fn=o{(l,Qi) = 1, and they largely simplify in the commutative chamber n — 00. 

There is a large freedom in choosing the value of k, and it is natural to ask if this choice 
has some physical interpretation. It was argued in |27] that this is indeed the case, and 
the choice of k is equivalent to the choice of open BPS chamber (open BPS chambers were 
introduced in section 12. ip . In particular, it turns out that the open generating parameter 
can be identified with matrix eigenvalues z^, and the open BPS generating function f lM|) 
in the open chamber labeled by k can be identified with matrix integrand 

Z^^;- = e-^^"("\ (78) 

Even though the overall factors in (177|) may involve closed moduli Qi, they do not involve 
open moduli z. In this sense the matrix integrand is well defined, and up to some simple 
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identification can be identified with open BPS generating function. This identification of 
parameters amounts to the shift z — )■ —zq^^'^ (to match earlier M-theory convention with 
half-integer powers of g, to integer powers of q in the fermionic formalism), as well as 
identification of Kahler parameters considered in M-theory derivation with parameters /ij 
introduced below. We also note that the BPS generating function in f lT^ is determined by 
the open topological string partition function associated to the external axis of the toric 
diagram, as in figure [lH As we will also see, the value of the above integral ( 177|) can be 
related to some more general Calabi-Yau geometry Y . 



/ 

/ 

/ 




Figure 14: Brane associated to the external leg of a toric diagram (of a conifold 
in this particular case). Closed string parameter is denoted by Q and open 
string parameter by z. 



5.2 Matrix models for the non-commutative chamber 

In this section we illustrate the relation between BPS counting and matrix models in 
case of the non-commutative chamber n = 0, and the choice of open chamber also k = 0. 
This corresponds to the insertion of the identity representation f l7^ exactly in between 
\Q±) states in (153|) . In this n = case no factor f^ in ([77D arises, and we obtain matrix 
models with potentials which can be expressed in terms of the following version of the 
theta function 

oo 

Q{z-q) = l[{l + zq^){l + q^+'/z). 

j=0 
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For a general geometry of the form shown in fig. |H1 with types of vertices given by ti, 
corresponding matrix models take form 

r r ^ 

Z= rff/JJe--^^^"^ = / dUlllle{U+^z^iq,---qiy,qy'^\ (79) 

''a a 1=0 

where integral is over unitary matrices of infinite size, N = oo. Special cases of this result 
include: 

• for C'^ the result fl79l) provides a matrix model representation of MacMahon function 
Z = M(l) = n^^il-'- ~ in terms of a matrix model of the form (17^ with the 
integrand 

oo 

e-^^(^) = + zq^){l + q^+'/z) = e{z; q). 

3=0 

• for the conifold we obtain a representation of the pyramid partition generating func- 
tion fll9p (with = 0) in terms of a matrix model with the integrand 



M (i + Q^g^)(i + ^) ©W^;^) 



• for C'^/Z^^^ we have tp = +1 for all p and we find matrix model representation of 
the BPS generating function in terms of a matrix model with the integrand 

fj^o {qi---qN)z 

N 

= ll'^{iqi---qi)z;q) 

1=0 

5.3 Matrix model for the conifold — analysis 

In this section we illustrate how matrix model techniques can be used in the context of 
models which arise for BPS counting. We focus on the conifold matrix model in arbitrary 
closed BPS chamber n, and fixed k = 0. In this case the result (1771) takes form (after the 
redefinition Q = —gig") 

oo 

Zn = M(g)2j](l-gg^y(l-g-V+T+"= (80) 

= fnM) J ^^nn (l + ,^,..n.l/g) (l + ,.g/,,) = /"(g'Q)^— 
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with 

with MacMahon function M{q) defined in (fTSl) . and with the following generalized MacMa- 
hon function 

i=l ^ ' 

In particular, in the non- commutative chamber = \^ a^fj in the commutative 

chamber f^^^^'^ = M{q) which represents topological string degree zero contributions. 
The result (IHOl) implies that the value of the matrix model integral (without the prefactor 
fn) is equal to 



^matrix 



fniq,Q) .^^ 



J\3 



where fi = q^/Q- 

Now we wish to analyze the matrix model Zmatnx- We parametrize the 't Hooft coupling 
and the chamber dependence respectively by 

T = g,N, T = ng,. (84) 

As our models correspond to f/(oo) matrices, ultimately we are interested in the limit 

T — !■ oo, Qs = const, Q = const, (85) 

for each fixed chamber (i.e. fixed n and therefore r). The non-commutative chamber 
corresponds to r = 0, while r — )■ oo represents the topological string chamber. 
Using the expansion of the quantum dilogarithm 



logn(l--^0=4ELi-4-)^> (86) 



i=l m=0 



and the redefinition of the unitary measure fl36|) we find, to the leading order in Qg, the 
following matrix model potential 

Vr = T \og{z) + Li^{-z) + Li2 ( - -) - Li2 ( - -) - Li2 ( - ^) , (87) 
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so that 

z 

Now we wish to solve the model f lHOj) in the small Qs limit. Firstly we need to find the 
resolvent uj{p), which can be done using the Migdal integral (137|) . and careful derivation 
is presented in [29j. As we expect one-cut solution of our model, from the Migdal integral 
we get an expression in terms of the end-points of this cut a and h. The normalization 
condition (138|1 imposes two constraints, for terms of order and p^^ in the resolvent, 
which take form 



VoTQ - y/b + Q ^1/2(t+T)/2 



Q^/^e^"+^^/' (89) 



V{a + Q)b- y/{b + Q)a ^ Qi/2^-{r+T)/2 .gQ^ 
^ia + Qe-)b- y/ib + Qe-] 



)a 

These constraints can be solved in the exact form, with result 



2 (1 - ^)(1 - /i6^) + (l-g) (l + /xe2 _ 2/i) 



where we introduced 



-i + e' ' " ' „^,z — —+ (91) 

+ (l-/ie2)2 

v/(l-g)(l-e2)(l-^)(l-g;xe2y 
(l-/ie2)2 



Substituting these end-points back to the formula for the resolvent we find 



. . 1 , ^ (1 + ^e^)p + (1 + Qe') T (1 - f^e')^ip - a){p - b) ^ 
-±(P) = ^ log ( 2e-^p + Q) ) • 



As a check, this result indeed satisfies the consistency condition fl39j) 

u+{p)+uj-[p) = — , 

with Vr given in fl87p . From the knowledge of the resolvent we can also determine eigenvalue 
density along the cut (HOj) 

1 / (I + ^e')p +1 + Qe^-{1- ^,e')^{p -a){p-b) \ 
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as well as the spectral curve. Writing x = pTuj{p) and p = e^, and after a few simple 
rescalings we find that the spectral curve takes form 



e-+y + e' + ey + Qi e^' + Q2 e^^ + = 0, (95) 

where 



2 l + f^Q 



;i + /.e2)(l + Q62)' 
1 + Qe' 

Qs = Q '^^'^ 



^l + e2Q)(l + /xg)' 




Figure 15: The spectral curve for the conifold matrix model / f^) in arbitrary 
closed BPS chamber coincides with the mirror curve of the closed topological 
vertex geometry, whose toric diagram is shown above. This geometry has three 
moduli, which encode conifold size Q, the closed BPS chamber via /x, and 
finite 't Hooft coupling via e~'^ . 



The above curve is given by a symmetric function of Q, fi = Q~^q" and = e~^. 
Apparently this is a mirror curve of the so-called closed topological vertex geometry, which 
is a Calabi-Yau manifold arising from a symmetric resolution of C'^/Z2 x Z2 orbifold, see 
fig. [151 This geometry has three moduli, i.e. the original Kahler moduli Q of the resolved 
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conifold, the chamber parameter n (encoded in ^) and the 't Hooft parameter T, which are 
all unified in a geometric way in our matrix model. Moreover, the fractional coefficients 
in the curve equation fl95|) encode the correct mirror map for the closed topological vertex 
geometry (and to the linear order, these coefficients are just Q, ^ and e~^). 

In the BPS counting problem we are interested in, as follows from the form of the 
identity operator 0761) . ultimately we need to consider matrices of infinite size. We also 
need to keep fixed Qs-, so we should consider the limit of T oo, or equivalently e — )■ 0. 
Up to a linear shift of x and the equation (195|) in this limit becomes 

/i + e"+^ + + (1 + Q/i) + Q = 0. (97) 

The manifold corresponding to this curve is the suspended pinch point (SPP) geometry, 
with Q and encoding fiat coordinates representing sizes of its two P^'s. Having found 
the SPP mirror curve, let us also make the following remarks. 

Firstly, we see that not only the spectral and mirror curves agree, but moreover the 
matrix integral (IHOj) reproduces (after the identiciation q = g^) the full topological string 
partition function of the SPP geometry at finite Qs, which is known to take form 



for Kahler parameters Q and /i. This confirms that our result makes sense, although 
this also means that the terms in lowest order in Qs in the potential reproduce the full Qs 
dependence of the partition function. It would be nice to prove rigorously that higher Qg 
corrections to the potential indeed do not affect the total partition function. This appears 
to be a very special feature of matrix models integrands which can be expressed - as is 
the case for f l83|) - in terms of inffnite products of the form Yli,{l — xq^). One proof of 
such statement in a related situation (although in addition taking advantage of a special 
phenomenon of the arctic circle) can be found in |52] . 

Secondly, it is natural to conjecture that the total partition function of the matrix 
model, for ffnite T, should reproduce (at least up to some MacMahon factor) the topological 
string partition function for the closed topological vertex which reads 



oo oo 



Qq^fil - iiq^'fil - t^q'^fil - Qiit'^q^f 



2nk\k 



Finally, we also note that in the limit Q, fi ^ our model reduces to the Chern-Simons 
matrix model discussed in [50] . Indeed, in this limit the potential (187|) reproduces 
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gaussian potential 



In this case the resolvent flMl) reduces to 



(p).-Q-o = ^ log ( 



and agreee^ with the resolvent of the old Chern-Simons matrix model found in [17t 150] , 
which is known to reproduce the Chern-Simons partition function 

_ M(g.) 

^Chern-Simons r / rp n ; I yy I 

M(e-^,gs) 

with 't Hooft coupling T identified with the Kahler parameter of the conifold. A manifes- 
tation of this result is seen directly in the product formula in (183|) . which is an expression 
for matrices of infinite size (i.e. T — )■ oo): for Q, /i — )■ the matrix model partition function 
reduces to the MacMahon function M{q), which is consistent with fl99l) in T — > oo limit. 
For finite T, also the spectral curve of the matrix model with the above gaussian potential 
reproduces the conifold mirror curve (HTl) of the size given by the 't Hooft coupling 

X + p + xp + x^e~^ = 0. 

Let us also mention that the MacMahon function arising in T — ?■ oo limit of this model, 
and in fact the entire Chern-Simons partition function, can be obtained - following the 
postulate of the remodeling conjecture [5,lj - from the topological recursion applied to the 
above curve [57j . It turns out that the topological recursion associates a square root of the 
MacMahon function to each patch of C^, in agreement with (1171) . and building the toric 
Calabi-Yau three-fold by glueing such patches is mirrored by the pant decomposition of 
the mirror curve. However, the behavior of such constant contributions to the partition 
function in matrix models is subtle, and in general may not agree with the topological string 
(or topological recursion) result, even if a spectral curve of a matrix model reproduces an 
appropriate mirror curve. For example, matrix models constructed in [55] reproduce 
correct mirror curves (in particular, for the conifold) and encode correct non-constant 
contributions to the partition function, however by construction do not involve any factor 
of MacMahon function. 



^Instead of introducing Tlogz term to the potential to get the standard Vandermonde determinant, 
the solution in [47, involves completing the square, which leads to a redefinition phere = P[|7]e^. Due to a 
different sign of gs we also need to identify 't Hooft couplings as There = ~i|47]- Taking this rescaling into 
account, our cut endpoints (191192^ with Q — fi = also coincide with those in [47]. 
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5.4 Matrix models and open BPS generating functions 

In this section we finally consider arbitrary closed and open BPS chamber, so that 
matrix models take most general form. Analyzing the case of C'^, conifold and C'^/Z2 
we illustrate the claim (ITHl) that open BPS generating functions can be identified with 
integrand of matrix models. Also for this reason our analysis is only on the level of these 
integrands, however it would also be interesting to understand the corresponding spectral 
curves. 

5.4.1 

We recall that the open topological string amplitude for a brane in is given by 
the quantum dilogarithm (IT^ . The condition for the central charge (II ip and the general 
formula (11^ imply that in the open chamber labeled by k the open BPS generating function 
reads 

oo oo 
^ope. ^ -Q^^ _ ^^.-1/2) -Q^i _ Z-\l-y'). (100) 

1=1 j>k 

In one extreme chamber A; = the generating function is equal (up to the overall g^/^"^) to 
a theta function, and so it is a modular form, as explained in section 12.11 On the other 
hand, for k — )■ oo, the generating function Zf^^oo reduces to the open topological string 
amplitude (fT3|) . 




Figure 16: Factorization of crystal which leads to open BPS generating 
functions. The size k encodes the open BPS chamber. 
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Now we present how the result fllOOp arises from the matrix model viewpoint. To start 
with, we again consider fermionic representation. Following results of section 14.4. in this 
case ^+(1) = r+(l)Q and to the geometry of C'^ we associate a state 

oo 

4=1 

and similarly for There is a single closed string chamber in which the generating func- 

tion Z = = M{q) is given by the MacMahon function. Following the prescription 

(I77|) we insert the operator I at the location k. This gives 

oo fe— 1 

Z = M(g) = (0|n^+(l)|Iin^+(l)l^-) = 

i=k j=0 

= f\q)Z^atrr.. (101) 

where the matrix integral is given by 



/oo oo 
a j=l i=k 



and 

~ ^ M(l) 



/'(«)=nnT3^ 



with the generalized MacMahon function M{q^, q) defined in (182|) . Matrix model integrand 
in Zmatrix indeed reproduces open BPS generating function (llOOp (up to a redefinition 
z — )■ —zq^l"^ and identifying q = g^) in a chamber labeled by k. 

5.4.2 Conifold 

Here we illustrate a relation between matrix models and open BPS generating functions 
for the conifold, related to a brane associated to the external leg of a toric diagram, as in 
fig- HH With appropriate choice of framing its amplitude reads 



-^1 



open -^(^5 Qs 



L{zQ,qsy 

This also reduces to the modular generating function in the non-commutative chamber 
n = k = 0. More generally, let us consider open BPS counting corresponding to the closed 
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Figure 17: Factorization of the conifold pyramid which leads to open BPS gen- 
erating functions. The size of the pyramid n represents the closed BPS chamber, 
while the size k encodes the open BPS chamber. 



chamber labeled by n, and open chamber labeled by k. From the condition (\TT\i . after the 
shift z — —zq^^"^, we get a general generating function of open BPS states 

2open,k _ \^open\2 _ TT + ^9'-;) ~^ ^ Qs^ ) (\[]'?\ 

- I top \ chamber " H Q ^ zOqiMl + Z-^O-^Q^+^+l-^) ' ^ ' 



zQqi){l + z-^Q-^qr^ 

This result arises from matrix model viewpoint as follows. We take advantage 
of the results of section I4.4.3[ where we determined A±{x) = T ±{x)QiT'^{x)Qo and 
W = r_(l)QirV(l)Qo- This leads to 

oo fc— 1 

Zn = (0| n^+(l)|I| n^+(l)l^l^-) = fn^rnatri.. (103) 
i=k j=0 

In terms of yU = — ^ = Q'^^q"' the matrix integral takes form 



■'matrix 



[l + z^q^){l + z-\'+^-') 



The integrand of this matrix model indeed agrees with (I102p M-theory (again identifying 
yU with Kahler parameter used in M-theory derivation, and setting q = qs)- In the limit 
n — )■ oo followed by — )■ we get the result for C'^ given in fllOip . On the other hand, for 
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both n, A; — )■ oo, the integrand reduces to the open topological string amplitude given by a 
ratio of two quantum dilogarithms. The prefactor above is found as 

fk ^ ,.( .2 M{^iq\q)M{Qq\q) . 

M(/x, g)M(g^ q)M{Q, q)M{^^Qq\ q) ^ ' ' 

5.4.3 CVZ2 

As another example we consider open BPS counting functions for resolved C^/Z2 sin- 
gularity. In this case the topological string partition function for a brane on the external 
leg is 

Z-* = L(z,g,)L(zQ,g,). 

This implies the following BPS generating functions in a closed chamber n and open 
chamber k (after z — )■ —zq^^'^ shift) 

00 

1=1 

On the other hand, using results of section |4.4.2[ i.e. A±{x) = r±{x)Qir±{x)QQ and 
W = r_(l)Qir+(l)(5o, and redefining /i = — = Q'^q"^, we obtain 
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Zk = (oin^+(i)iiin^+wi^i^-)=/'^" 

i=k j=0 

-.fc+jr — 1 „n+k+j~l 



^uHHii + + W)(i + )(i + )• (104) 

a j=l 

The matrix integrand indeed agrees with the M-theory result (when written in terms of 
the argument /i) above. The prefactor above reads 

fk ^ j^(.2 M{fi,q)M{Q,q) _ ■ 

M(/xg^ g)M(g^ q)M{iiQq\ q)M{Qq\ q) ^ ' 

6 Refined crystals and matrix models 

In the last section we turn our attention to so-called refined BPS amplitudes, and 
explain how to incorporate the effect of such refinement in the fermionic formalism and 
matrix models, following [21]. To start with, we recall that there are various definitions of 
refinements, which arise in the context of BPS counting or topological string theory. Here 
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we focus on closed BPS states and consider the following characterization. We introduce 
an additional parameter y on which multiplicities of D6-D2-D0 states f2 in the original 
definition of the generating function ([1]) may depend 



For fixed DO-brane and D2-brane charges a and 7, and a choice of closed BPS chamber n, 
refined degeneracies are defined as 



where 'Ha,-y{n) denotes a space of BPS states with given charges a, 7 and asymptotic values 
of moduli corresponding to a chamber n, and J3 represents a generator of the spatial rota- 
tion group. For y = 1 these degeneracies reduce to those in ([1]). These refined degeneracies 
are interesting invariants if the underlying Calabi-Yau space does not posses complex struc- 
ture deformations - and this is indeed the case for non-compact, toric manifolds we are 
interested in. In jlOj it was argued that these invariants agree with motivic Donaldson- 
Thomas invariants of |2], and in the case of the resolved conifold the corresponding BPS 
generating functions were derived using the refined wall-crossing formula, and encoded in 
a refined crystal model. From mathematical viewpoint this setup was generalized to the 
whole class of toric manifolds without compact four-cycles in [13], and shown therein to 
agree, in the commutative chamber, with refined topological vertex computations. The re- 
fined topological vertex was introduced in [12], see also |58 | [59 | 160 ] . For other formulations 
of refinement see |Tl] |6T| 162] . 

Our aim in this section is to construct refined crystal and matrix models, which would 
encode refined BPS generating functions, and in particular (in the commutative cham- 
ber) refined topological string amplitudes. We note that an additional motivation to find 
such models arises from the AGT conjecture [63] and the results of |64] . which state that 
partition functions of Seiberg-Witten theories in the ^-background (which are related to 
topological strings by geometric engineering) are reproduced by matrix models with /3- 
deformed measure (i.e. with Vandermonde determinant raised to the power /3). Explicit 
construction of one class of such /3-deformed models, however only to the leading order, was 
given in ^55j; some other works analyzing five-dimensional beta-deformed models include 
[65| [66] . On the other hand, explicit computations for simpler /3-deformed model with 
gaussian potential [67t [68] . revealed that it does not reproduce refined topological string 



Z:^f{q,,Q) = J2K:^in;y)qfQ\ 




(105) 
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amplitude for the conifold (even though the unrefined topological string partition func- 
tion is properly reproduced when /3 = 1, see |17t I5U]). Nonetheless, the question whether 
there exist matrix models which reproduce such refined amplitudes remained valid. As we 
show below (following |3T]) such models can indeed be constructed by appropriate defor- 
mation of the matrix model potential (rather than the measure). We note that recently 
another class of matrix models (with different than above deformation of the measure) was 
proposed [69], which also reproduce refined generating functions. 

Let us also note that in this section we consider the same set of walls as in the unrefined 
case. More general walls, along which only refined BPS states decay, may also exist [lOj. 
They are called invisible walls and they do not arise in our analysis. 

In this section we use the following refined notation. The string coupling Qs, related to 
the DO-brane charge as = e~^* in the unrefined case, is replaced by two parameters 



or equivalently (3 = —^,€162 = —g^. We also often use the exponentiated parameters 

ti = e-'\ t2 = e'\ (106) 

and introduce 

gsB = 61 + 62 = gs{\/l3 - ■ 

The variable y in fllOSp is related to ti and t2 as y = ti/q^ = qs/t2, so that y"^ = ti/t2 = ■ 
In this notation the unrefined situation y = 1 corresponds to (3 = 1, for which ei = —62 = gs 
and ti = t2 = Qs and B = 0. 

Let us present now refined BPS generating functions for some Calabi-Yau spaces: 

• For C'^ we get the refined MacMahon function |12] 

00 ^ 

= M.e/(ti,t2) = n 

k,l=0 ^ ^1 ^2 

In this case there is no Kahler parameter, and therefore there are no interesting 
wall-crossing phenomenajf] 



^In fact one can consider more general family of refinements parametrized by S, such that Ms{ti, 12) ~ 
nfe°/=o {^^^1 ^ ^2 ^ ) • tl^i*^ paper we fix the value 6=1 (note that in [fOj another choice 6 = 
was made). 
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• For the resolved conifold refined generating functions were computed in [TU] using 
refined wall-crossing formulas. In the closed BPS chamber labeled by n — 1 these 
generating functions read 

oo 

^conifoi, ^ M^,f{t„hf[ n (1 - Qti^%)) ( n (1 - Q-'t'A))- (108) 

k,l=0 k>l, l>0, k+l>n 

In the commutative chamber — )■ oo the terms in the last bracket decouple and 
the BPS generating function agrees (up to the refined MacMahon factor) with the 
refined topological string amplitude 

oo 
k,l=0 

On the other hand, in the non-commutative chamber n = the refined generating 
function is given by the modulus square of the refined topological string amplitude. 

• For a resolution of C^/Z2 singularity there is also a discrete set of chambers 
parametrized by an integer n. The corresponding BPS generating functions read 

oo 

= M..,(ti,t2)^( n {l-Qt'^^%r') ( n {i-Q-'t^y). (m) 

k,l=0 k>l, l>0, k+l>n 

• Generating functions for an arbitrary toric geometry, in for the non- commutative 
chamber, are given (as in the unrefined case) by the modulus square of the (instanton 
part of the) refined topological string amplitude 

Zo'' = \zrd? = zr4mzroi{Q-'). (no) 

The (instanton part of the) refined topological string amplitude is given by [T2l [58] 

oo 

Zi:l{Q,)=M,,j{t,M)'^X{ n [^-{Q^Q^+l■■■Q^-l)ttXy"^'■ (111) 

fc,Z=0 l<i<j<N+l 

6.1 Refining free fermion representation 

In the non-refined case to a geometry consisting of P^'s we associated in section H] a 
crystal which can be sliced into layers in + 1 colors, denoted go? O'l? • • • , Qn- In that 
case parameters qi, . . . ,qN encode Kahler parameters of the geometry Qi, . . . ,Qn, while 
the product Hilo^* mapped to (possibly inverse of) = e~^^ In the refined case the 
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assignment of colors must take into account a refinement of a single parameter into 
ti and ^2 introduced in (I106p . In particular, in the non-commutative chamber gj^o are 
mapped (up to a sign, as in the non-refined case) to Qi, however we will have to replace 
go by two refined colors q^^^ or ql^\ so that ti = Qq^i ■ ■ ■ Qn, for i = 1,2. The simplest case 
of C'^ refined plane partitions (discussed also in jT2]) is shown in fig. [TSi In what follows 
we will discuss assignment of colors for other manifolds. 




Figure 18: Refined plane partitions which count D6-D0 hound states in <C^ . In 
each slice halls, which intersect a dashed or solid line, have respectively weight 
ti or t2- The resulting generating function is the refined MacMahon function 

Mref{tiM)- 



Now we wish to follow the idea of section |H Firstly, we wish to construct refined states 
l^±'^) whose correlators would reproduce refined BPS amplitudes in the non-commutative 
chamber 

zi"^ = (112) 

Secondly, we wish to construct refined wall-crossing operators , such that the BPS 
generating function in n'th chamber can be written as 

Z;^/ = (113) 

In section [6.1.11 below we construct states \^±'^) for arbitrary manifold in a class of our 
interest. In section [6.1.21 we construct states \^±'^) and wall-crossing operators W^^^ for 
all chambers of the resolved conifold and a resolution of C'^/Z2 singularity. 
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6.1.1 Arbitrary geometry — non-commutative chamber 

Here we construct fermionic states which aUow to write the BPS generating 

functions in the non-commutative chamber as in flll2p . As in the non-refined case, the 
states \^±'^) are constructed from an interlacing series of vertex operators and weight 
operators. The refinement does not modify the three-dimensional shape of the correspond- 
ing crystal, therefore the assignment of vertex operators is the same as in the non-refined 
case (H3|) . as explained in section 14.2. II However, this is assignment of colors, encoded in 
the weight operators, which is modified in the refined case. Let us introduce N operators 
Qi representing colors gj, for i = 1, . . . , A, and in addition two other colors ^q^^ and ql^\ 
which are eigenvalues of qI^^ and Q^^ ■ Operators Qi, . . . ,Qn, similarly as in section |4. 2. ![ 
are assigned to P^'s in the toric diagram, and we introduce 

Q^'^ =Qil^Qi---QN, U = qfq^---qN, for z = 1,2. (114) 

Now we define refined version of A± operators 

A_{x) = r:nx)gir!?(x)Q2 

Commuting all Qj's to the left or right it is convenient to use 

qi qiq2 qiq2 ■ ■ ■ qN 

A_{x) = A_{x) {Q^'Y' = rL^(x)r!?(xgi)r!!(a;gig2) ■ ■■Tj^'{xq,q2 . ..qN), 
and when the argument of these operators is a; = 1 we often use a simplified notation 

A± = Z±(l), A± = A±(l). 

Finally we associate to a given toric manifold two (refined) states 

= (0|...A+(l)Z+(l)A+(l) = (0|...A+(t?)A+(ti)A+(l), 

in:!^) = Z_(i)A_(i)Z_(i)...|o) = A_(i)A_(t2)A_(t^)...|o), 

where |0) is the fermionic Fock vacuum. 

Our claim now is that the refined BPS generating function can be written as 

Zo^^^ = {n^_^^\QL'f) ^ Z,,,{Q,)Z,^{QT^), (115) 



■■■TI-{x)QmTT\^M\ 
■■■r-{x)QNTj^\xM^. 
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with ZtopiQi) given in fl23p . and if one identifies gj parameters which enter a definition of 
\^^±^) and string parameters Qi = e~"^' (for i = 1, . . . , N) as follows 

and with refined parameters ti,2 identified as in flll4p . 

To prove flllSp for general geometry, we first note that commuting operators A^{x) 
with A^{y) 

A+{x)A4y) = A4y)A+{x) C{x,y), 

gives rise to a factor 

Now we write the states in terms of A± operators, and commute T± within each 

pair of A+ and A_ separately 

oo oo oo 

Zo'^ = = (oi(n^+(n)) (n^^(^2))io) = n ^'(n,^^- 

i=0 j=0 i,j=0 

This last product reproduces modulus square (11151) of the refined topological string parti- 
tion function (lllip . and therefore proves the claim (I112p . Moreover, for the special /3 = 1, 
we automatically obtain the proof of the analogous statement (13^ in the non-refined case 
from section 14. 2. 21 

6.1.2 Refined conifold and CVZa 

We can now extend the fermionic representation to non-trivial chambers, for simplicity 
restricting our considerations to the case of a conifold and a resolution of C^/Z2 singularity, 
which both involve just one Kahler parameter Qi = Q. Our task amounts to determining 

re/ 

appropriate wall-crossing operators, denoted so that in the chamber labeled by n — 1 

the BPS generating function can be written as 

Z:-{ = (^^+ Wn-ll^-')- (116) 

In these both cases the toric diagram has two vertices, the first one of type ri = 1 and the 
second one denoted now r = T2, and r = =pl respectively for the conifold and C^/Z2. A 
crystal associated to the expression (I116p has n stones in the top row and can be sliced into 
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Figure 19: Refined pyramid crystal for the conifold, in the chamber correspond- 
ing to n stones in the top row. Along each slice (as indicated by broken or 
solid lines) all stones have the same color, assigned as follows. On the left side 
(along broken lines), each light (yellow) and dark (red) slice has color denoted 
Qq and qi respectively. Moving to the right, in the intermediate region (along 
solid lines), a color of each new light or dark slice is modified by respectively 
q^^ factor (with respect to the previous light or dark slice). On the right side 
(again along broken lines), each light or dark slice has again the same color, 
respectively qoq^^" or qiq^'^^^^\ The assignment of colors in the intermediate 
region (along solid lines) interpolates between constant assignments on the left 
and right side of the pyramid. 



interlacing single-colored layers. The assignment of colors is analogous as in the pyramid 
model discussed in [101 [13] • The pyramid crystal for the conifold and C'^/Z2 are shown in 
fig. [I9]and|20l 

The assignment of colors is determined as follows. All stones on one side of the crystal 
are encoded in 

{n^^^l = (o| . . . (r+(i)Qir;(i)Qo) (r+(i)Qir;(i)Qo) • 

The Kahler parameter Q, as well as the parameter ti, arise as 

qi = rQtl'"-, qQ = r-^, so that qoqi = ti. 
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Now the crystal with n — 1 additional stones in the top row arises from an insertion of 
the operator 

w":!, = (r_(i)gir;(i)goP)(r-(i)Qi?r;(i)gof^)... 
. . . (r_(i)Qig(^^r;(i)go^^^). 

This operator is made of n — 1 terms of the form (t ^{l)Qiq^^T'^{l)QQq^^^+^^^^ for i = 
0, . . . , n — 2, where in each subsequent dark or light slice we insert one additional operator 
q^^. This additional operator changes the weight of each stone in this slice by q"^^ = 
{ti/t2)'^^ (with respect to the previous slice of the same light or dark color). 

Finally, all stones on the right side of the crystal have again the same light or dark 
color, so that the corresponding state is 

= (r_(i)gig(^^r:(i)gof^) (r4i)gig(^^r:(i)gof^) . . . |o). 

We see that varying weights in the middle range (along solid lines in fig. [19] and |20|) 
interpolate between fixed weights of light and dark stones on two external sides of a crystal. 




Figure 20: Refined pyramid crystal for the resolution o/C^/Z2 singularity, in 
the chamber corresponding to n stones in the top row, as seen from the bottom 
(i.e. a negative direction of z- axis). Even though the three-dimensional shape 
of the crystal is different than in the conifold case, the assignment of colors is 
the same, see fig. [73 
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We can now commute away all weight operators in the above expressions, using com- 
mutation relations from section 13.21 This results in 

oo n—2 oo 

ZnJ, = (o| ( n r+(t?)r;(t^M)) ( H r-(4)r;(gr'ir)) ( H r-{tr'^')rL{tQ4)) |o). 

k=l i=0 k=0 

(117) 

To check that this is a correct representation we commute all vertex operators, and find 

oo oo 

ZZ{ = Mref{h,hr n i^-Qt'A)-^ n il-Q-'t%)-\ (118) 

fc=l,«=0 k>l,l>0,k+l>n 

where r = =1=1 respectively for the conifold and C^/Z2. This result reproduces fllOSp and 
fllOOp . which confirms that the fermionic representation we started with is correct. 

6.2 Refined matrix models 

In the refined case one can associate matrix models to refined generating functions 
in the same way as described in section 15. H i.e. by inserting the representation (I76|) of 
the identity operator into fermionic representation 01151) or flll6p . This does not change 
a unitary character of the matrix model, which is a consequence of the representation 
(176|) . However, due to more subtle weight assignments, this is matrix potential which gets 
deformed by /3-dependent factors. In general we will therefore obtain matrix models of the 
following form 

Z-^ = fn [vulle''^''^'-'^\ (119) 
k 

where for convenience we introduced a factor ^/P in front of the potential V{z; /3). We will 
consider a few examples below. 

6.2.1 Non-commutative chamber 

For arbitrary geometry, in the non- commutative chamber, refined matrix model inte- 
grand can be expressed in terms of the refined theta function 

oo 

e{z;h,t2) = l[{l + zt{-''){l + ti/z). 

Repeating the computation described in section 15. Ij however starting with the refined 
representation flll5p . in the non- commutative chamber for general geometry we find the 
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matrix model 

i.e. we identify e~~a^^'^^'^^ = Ililo ® ' ' ' 5'/) ""^^ ^i; ^2)^'^^- The product over / runs 

over all vertices and we identify Kahler parameters Qp with weights qp via qp = {TpTp^i)Qp. 

6.2.2 Refined matrix model 

We obtain a refined matrix model for as the special case of f ll20p . For the refined 
the BPS generating function is a refined MacMahon function Z"^^^ = Mref{ti, ^2) introduced 
in (I107p . and the corresponding matrix integrand takes form of a a refined theta function 

e-f ^ Ylil + ^ti+^)(l + i/z) = eiz; t„ h). (121) 

j=0 

Using the asymptotics (186|) we find the leading order expansion of the potential 

g-^V{.;/3) _ [_i(log.)2-(l-/3-i)Li2(-^)+0{s.,/3)] _ ^^22) 

The quadratic term in the potential is the same as in the non-refined case. The term 
involving Li2(— 2;), as well as all higher order terms 0{gs,P), arise as deformations which 
vanish for /3 = 1. Therefore, for /3 = 1, we obtain a Chern-Simons matrix model which 
indeed gives rise to MacMahon function in — )■ 00 limit, as we explained in section 
15.31 For arbitrary /3, the resolvent uj{p) can also be found using the Migdal integral (137|) . 
In principle one could repeat the computation described in section 15. 3[ however this is 
technically more involved. Nonetheless, this would lead to /3-deformed end-points of the 
cut (IHT]) and and in consequence to the /3-deformed spectral curve. This curve would 
be some /3-deformation of the mirror curve given in (HTl) . It is still an interesting question 
to find this curve in the exact form and analyze its properties. 

6.2.3 Refined conifold matrix model 

Finally we find matrix models for the refined conifold. Starting with the represen- 
tation (I117p . inserting the identity representation (176|) and following standard by now 
computations, we find the following matrix model for the conifold in the n'th chamber 
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(corresponding to a pyramid with (n + 1) stones on top) 

oo oo 

= M,e/(ti,t2)' n n i^-Q^'tit 

k=l,l=0 k>l,l>0,k+l>n+l 

+ (l+ti/Zk) 



with the prefactor given by 



Ti OO 71 OO 

/n(g,Q)=(nnrzW)(n n a-n^^^/Q) 

j=l fc=0 ^ 2 j=l j=n+l-i 

In the hmit of the commutative chamber, n — t- oo, we get = Mref{ti,t2). Therefore in 
the commutative chamber we get a matrix model representation of the refined topological 
string conifold amplitude 

oo 
k,l=0 

i^ + tiQ/zk) 

In this case the lowest order potential is a modification of the potential fll22p by a 
Q-dependent dilogarithm term 

Viz; /3) = -^(log^)2 - (1 - p-')U,i-z) - Li2(-g/^) + Oig^, /?). (123) 

This is quite an interesting result - as we already explained above, it has been postulated 
for some time that the refined topological string amplitude for the conifold should have 
matrix model representation, however it was not clear how to derive it. Here we find an 
explicit matrix model representation of this amplitude. The corresponding spectral curve 
would again be a /3-deformation of the conifold mirror curve fl4ip . It would be interesting 
to compare it with other notions of deformed, or quantum mirror curves in literature. We 
also note that in the limit Q — )■ the above topological string partition function becomes 
just the refined MacMahon function, and the matrix integral consistently reproduces 
result (fm]) . 
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